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Abstract. We use the Gribov-Zwanziger Lagrangian in QCD to evaluate the one loop correction 
^ ■ to the gluon propagator as a function of the Gribov volume and verify that the propagator 

vanishes in the infrared limit. Using the corresponding correction to the Faddeev-Popov ghost 
propagator we construct the renormalization group invariant coupling constant, a^f(p 2 ), from 
the gluon and ghost form factors and verify, using the Gribov mass gap condition, that it 
freezes out at zero momentum to a non-zero value. This is qualitatively consistent with other 
approaches. We also show that there is an enhancement of the propagator of one of the Zwanziger 
ghosts at two loops similar to that which occurs for the Faddeev-Popov ghosts. From the exact 
evaluation of the form factors we examine power corrections for the gluon propagator and the 
effective coupling. We find that both have the same qualitative behaviour in that the leading 
power correction is 0(l/p 2 ) and not 0(l/(p 2 ) 2 ). 
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1 Introduction. 



In the late seventies, Gribov pointed out that fixing a gauge in a non-abelian gauge theory was 
a non-trivial exercise in comparison with an abelian theory, PQ. In particular the gauge field 
was not fixed uniquely since for some gauge fields one could construct at least one other copy 
which could be obtained from the first by a gauge transformation. To understand this problem 
and the relation to confinement, Gribov restricted the region of integration in the path integral 
defining the non-abelian field theory to the region bounded by the first zero of the Faddeev- Popov 
operator, pQ. The boundary of this region, known as the Gribov horizon, defines the Gribov 
volume and is parametrised by the quantity 7, of mass dimension one, known as the Gribov 
mass. Though inside the Gribov horizon one can still have gauge copies, [2 El II]- However, 
within this horizon is the fundamental modular region where each gauge field is uniquely defined, 
without copies, upon gauge fixing. Moreover, it had been suggested, p], that any additional 
gauge copies inside the Gribov horizon do not affect the vacuum expectation values of any 
operator. Subsequent zeroes of the Faddeev-Popov operator beyond the first, define further 
regions of configuration space but it has been shown, pQ, that they are equivalent to the first 
Gribov region which includes the origin. If one accepts that quantum chromodynamics (QCD), 
which is clearly a non-abelian gauge theory, has to be restricted to the volume defined by the 
Gribov horizon then Gribov argued that this had a profound effect on the infrared structure of 
the field theory, [HI El 13 El El EH E] ■ For instance, in the Landau gauge Gribov showed that the 
ghost propagator diverged as jjsjz as p 2 —* rather than the usual ^ behaviour which occurs 
in the ultraviolet (and perturbative) region, pQ. This property, known as ghost enhancement, 
followed directly as a result of the Gribov mass satisfying a gap equation which was determined 
at one loop explicitly, 1 . A comprehensive intoduction to the original Gribov article is provided 
in [TU. 

In recent years, given the advance in non-perturbative techniques such as Dyson-Schwinger 
equations (DSE), the behaviour of the ghost propagator has been studied in the Landau gauge 
in the infrared and is in agreement with Gribov's enhancement observation, |131 1141 1151 116j . 
Though the actual power law exponent differs from that of Gribov, |131 1141 ITdI I16j . Further, 
such studies have opened up the possibility of examining other infrared features. For instance, 
it is believed the gluon propagator does not diverge at zero momentum, as its ultraviolet form 
would suggest, but either tends to zero or a finite value. Whilst Gribov demonstrated the gluon 
propagator vanished at the tree or classical level in the Landau gauge, , there appears to be no 
definitive agreement which of the latter scenarios persist in the quantum theory. For instance, a 
recent lattice study argues that the data imply a finite non-zero value for the gluon propagator 
at zero momentum, |17j . Whilst stated in these general terms it is worth bearing in mind that 
there are technical considerations such as gauge variance which may result in different pictures, 
when viewed in different calculational approaches and renormalization schemes. Though we 
note that throughout, our remarks will always be concerned with the Landau gauge. Equally 
controversial is the behaviour of the strong coupling constant as(p 2 )- Phenomenologically it is 
believed that as(p 2 ) freezes out at zero momentum to a finite non-zero value, [T8l H9 l l20 | l2 H 
1221 1231 1241 1251 1261 1271 128j . Though as we will illustrate later what its precise value is appears 
to depend on the method used to extract it from experimental data as well as how one chooses 
to define it. More theoretical predictions based on properties of the underlying quantum field 
theory itself such as analyticity arguments or DSE, rather than the fits to experimental data, 
favour a value higher than those from data, |29| I3()| 13 1[ 1321 133| . Again these comments need to 
be tempered with the caveat that theoretical predictions are made as a consequence of certain 
approximations and in renormalization schemes which are mass dependent, which imply a gauge 
dependent coupling constant. However, as the quantity measured in DSE and lattice studies is 
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a renormalization group invariant in the Landau gauge, which is measured at zero momentum, 
the hope is that such approximations do not significantly affect the estimates. Though some 
lattice and DSE studies do not see a finite freezing but instead observe an effective coupling 
which vanishes at zero momentum, [341 1351 136j . Whether this is due to finite volume effects as 
argued in [HH] by DSE studies, but subsequently rebutted in [S3 is still not clear. Moreover, it 
would seem that one of the issues in this technical debate centres on the validity and application 
of the underlying Slavnov- Taylor identities to the fully non-perturbative regime which governs 
the infrared behaviour. Though it is worth noting that an effective coupling can be defined from 
any of the vertices of the QCD Lagrangian and the infrared behaviour of each does not have to 
be the same. 

In the context of the Gribov approach to the infrared structure the advance made by 
Zwanziger, [7|l8lllf)|lll|. in reformulating the Gribov path integral and its inherent non-locality 
deriving from the limitation to the Gribov region is significant. In particular Zwanziger con- 
structed a localized Lagrangian which involved extra ghost fields (which we will refer to as 
Zwanziger ghosts in contradistinction to the Faddeev-Popov ghosts), which when eliminated by 
their equations of motion reproduces the Gribov formulation of the problem. As of nearly equal 
significance is the fact that this Lagrangian (which we will refer to as the Gribov- Zwanziger 
Lagrangian) is renormalizable, jlUl 1381 139j . and explicitly contains the Gribov mass 7. This 
former property has been analyzed through the algebraic renormalization machinery and the 
Slavnov- Taylor identities have been established, [SHUSHl- Interestingly despite the presence of 
the extra fields and a mass, no new independent renormalization constants are required in the 
Landau gauge to render the Gribov-Zwanziger Lagrangian finite. Moreover, as is evident in 
QCD when massive quarks are present the MS renormalization of the coupling constant, gluon 
and Faddeev-Popov ghost fields remain unaffected by a non-zero 7, [SHI EH]- Crucially the gen- 
eral formalism implies that the Gribov-Zwanziger Lagrangian can be used for performing loop 
calculations. Indeed Zwanziger, [S], reproduced the one loop gap equation of [Tj by evaluating 
the vacuum expectation value of the Gribov mass operator which is equivalent to the implemen- 
tation of the Gribov horizon condition. More recently this has been examined at two loops in 
[10] where the two loop correction to Gribov's gap equation has been determined in the Landau 
gauge with massless quarks in the MS scheme. As a check that the gap equation is consistent, 
the two loop correction to the Faddeev-Popov ghost propagator was computed and shown that 
the ghost propagator enhancement is preserved as a behaviour precisely as a result of the 
two loop gap equation, jlOj . 

Thus armed with a renormalizable localized Lagrangian, which at two loops appears to 
reproduce the observed behaviour noted by Gribov, we can attack the problem of the behaviour 
of other quantities in the infrared which we have discussed above. This is the aim of this article 
where we will compute the one loop correction to the gluon, Gribov-Zwanziger and Faddeev- 
Popov ghost propagators exactly. As will be evident this is not the straightforward exercise it 
would normally be in the absence of the Gribov-Zwanziger ghosts or the Gribov mass. However, 
we will be able to study the p 2 — > behaviour of the gluon propagators and see if it vanishes at 
the one loop level or not. Equally worth studying is the effective coupling constant freezing at 
zero momentum. Since lattice and DSE studies examine an effective coupling constant which is 
renormalization group invariant but which depends only on the gluon and Faddeev-Popov ghost 
form factors, we will be able to construct this quantity and then take the p 2 — > limit. We 
believe this will be constructive for reasons beyond the qualitative calculation we will perform. 
Specifically, previous studies have primarily been numerical and deduced either from experiment 
or by solving QCD with lattice or DSE techniques. The actual mechanism and which of the 
quantum fields truly drive the freezing have not really been studied. Having a Lagrangian 
where the calculations can be constructed explicitly will, we believe, be important in furthering 
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that understanding. Given that we evaluate the propagators as exact functions of 7 and the 
momentum, we can also study the corrections to the effective coupling in the alternative limit 
of 7 2 — > 0. The motivation for this is to examine whether one can not only produce power 
corrections but also to see whether such corrections are 0(1/ p 2 ) as suggested by the lattice 
computation of jJsHH] or 0(1/ (p 2 ) 2 ) as would be expected on the grounds of a condensation of 
the gauge invariant operator G^ v G a ^ v . Finally, in referring to the work we present as qualitative, 
it is worth briefly clarifying at the outset that what we mean by this is that performing one loop 
calculations, where behaviour in the infrared appears to be consistent with non-perturbative 
highly intensive DSE and lattice studies, our analysis should only be regarded as comparative. 

The paper is organised as follows. In section 2 we recall the main structure of the Gribov- 
Zwanziger Lagrangian and discuss the propagators and gap equation. The formal background 
to the computation of the one loop corrections to the 2-point functions of the Gribov-Zwanziger 
Lagrangian is introduced in section 3 prior to giving the details of the explicit computation in 
section 4. Section 5 is devoted to the problem of the freezing of the effective coupling constant 
in the p 2 — > limit whilst we discuss the structure of the propagators and effective coupling 
constant in the other limit r y 2 — > in section 6. Finally, our conclusions are given in section 7 
and an appendix collects the full expressions for the exact 2-point functions. 



2 Gribov-Zwanziger Lagrangian. 

We begin by defining the localized Lagrangian derived by Zwanziger, jH], which replaces the 
non-local formulation of the gauge fixing condition of Gribov, PP, concentrating on the Landau 
gauge. In addition to the usual Faddeev-Popov ghosts, c a and c a , of the canonical gauge fixing 
construction, the Gribov-Zwanziger Lagrangian involves additional Zwanziger ghosts. These 
fields, denoted by and {uf,wf}, are spin-1 and carry two colour indices. The latter 

set are anticommuting Grassmann fields whilst the former are commuting. Thus the full Gribov- 
Zwanziger Lagrangian is, jS], 

L GZ = L QCD + $*Hff> (p^y* _ u abfl d v {D I/ to^) ab 

- gf^U^^cf^ - jf^Aa^bc + f abc A a,^ + ^^(2.1) 

with 

L QCD = \ G a G a»v _ (QtiAlf - c^d^D,c a + ift 1 'flip* 1 (2.2) 

where is the gluon, ijj is the (massless) quark and the indices lie in the ranges I < a < 
1 < i < Np and 1 < I < Nt where Np and Na are the respective dimensions of the fundamental 
and adjoint representations and Nf is the number of quark flavours. The covariant derivatives 
are given by 

D^c a = - gf abc Ay , D^ u = d^ 1 + igT a A a ^ u 

{D^ v ) ab = d^f - gr A A^ db . (2.3) 

For completeness, we have included the covariant gauge fixing parameter a which is in principle 
required in the derivation of the gluon propagators. The Gribov mass, 7, is present in the 2-point 
term mixing between the gluon and the Zwanziger ghosts. In various articles by other authors 
the factor of -7= is absent. However, as we will observe later this choice is essential to retaining 
a gluon propagator which is the same form as that originally derived by Gribov, pQ. Though, 
we note now that changing the numerical coefficient of 7 does not alter the overall physics but 
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merely rescales the value of the gluon mass as a function of 7. The presence of the mixed 2-point 
term is problematic from the point of view of developing Feynman rules since one would have 
a mixed gluon Zwanziger ghost propagator. One way to circumvent this would be to redefine 
the and § ab fields in such a way that the 2-point term is absent. It transpires that this is 
cumbersome to implement and in earlier work, [301) an approach has been found which allows 
one to systematically evaluate Feynman diagrams with a mixed propagator and to derive the 
two loop correction to the Gribov mass gap. 

To derive the Feynman rules from 1)2.1 J) for the vertices of the interaction Lagrangian is 
straightforward. However, the main difficulty is to determine the propagator for the gluon in 
the presence of the mixing term. Although this was described in [301 ^ ^ s instructive to develop 
the derivation since it serves as the basis for deducing the one loop correction to the gluon form 
factor. First, we note that in general in deriving the propagators from the quadratic part of a 
Lagrangian the fields are first transformed to momentum space. Therefore, concentrating only 
on the gluon and Zwanziger ghost part of ()2.1|) we have 



L 



GZ 
quad 



2^("P) 



1 

1 — 

a 



Kip) 



l_ r tc A a { _ p)4)b c, {p) 



+ —f 



aba 



H-pR(p) 



v 



2 **"(-p)<(p) • 



(2.4) 



If 7 was zero then the inversion of the momentum space operator associated with the terms 
quadratic in the fields proceeds in a textbook manner. However, for 7 non-zero we must first 
write L^ & ^ in matrix form with respect to the basis The inclusion of the factor 

■A= into the basis is crucial since L^ & ^ clearly now becomes 



(l-i) 

"7 f cab ^ u 



"7 2 f acd r] t11 ' 



-P 



a?Ap) 



Inverting the matrix using the unit matrix 



(2.5) 



then the a dependent propagators are 



be I 



icd 1 



p)> = 


- 5 ab p 2 


p)) = 


jabc^l 


p)) = 


5 ac 5 bd 


r.2 




[(p 2 ) 2 + Cai a ] 
pAp) 



+ 



pAp) 



+ 



[( P 2 ) 2 + aC A 7 4 ] 



P 



[(p 2 ) 2 + aC A7 4 ] 
aL^(p) 



[(p 2 ? + Cai a ] 



+ 



where 



PfiPv 



Lfivip) 



[(p 2 ) 2 + a^ 7 4 ] 

PliPv 



(2.6) 



(2.7) 



(2. 



p" p* 
are the respective transverse and longitudinal projectors and we have been careful to restore 
the basis vector j^^i an d its transpose prior to deducing ()2.7)) which is the origin of the 
location of the factor of — ?= in the mixed propagator. Further, without the factor of —7= in the 
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mixed term of (|2.1j) one would not only have \pl^ 2 in the off-diagonal terms of the matrix of 
(|2.5I) but also one would have a common factors of [(p 2 ) 2 + 2C A r y 4 ] and [(p 2 ) 2 + 2aC J 47 4 ] i n the 
propagators. In the Landau gauge calculation of (Tj this additional numerical factor of 2 was 
absent and we choose to be consistent with [Q as our convention. Moreover, as we will later 
have to write [(p 2 ) 2 + C A j 4 ] as the product of the more conventional factors \p 2 + i\/C A j 2 ] and 



[p 2 



1 Ca7 2 ] the choice in (|2.1[) therefore also turns out to be the most convenient. Next we 



note that in the Landau gauge the final propagators for the gluon and Zwanziger ghosts are, [H] , 



(A a Jp)Al(-p)) 



§ ab p 2 



2\2 



(A»e(-P)> 



«(p)#"(-p)> 



cd i 



V2[(p 2 ) 2 + C A ^] 

§ ac § bd 



-pAp) 



(2.9) 



r ■ p 2 [( P 2 ) 2 + c A ^y 

Whilst this may be an unusual way of considering the quadratic terms we note that it is primarily 
forced, aside from the original mixing, by the fact that the gluon is a real field whereas the 
Zwanziger ghosts are complex. Finally, for completeness we note that the propagators of the 
remaining fields are 



(c a (p)c?(-p)) 



«(p)^(-p)> 



6 ab 
P 2 

p2 



{A{PW J {-P)) = iS ij 5 IJ ^ 



(2.10) 



Having summarized the derivation of the mixed propagators we note that the § ab fields are 
the ones which implement the Gribov horizon condition, jS]. To recap from (|2.1I) the equation 

n g ives 

1 



of motion of 



bf 



7 rabc 

^2 J ~d v D v "» 



-A c , 



whence one has 



f abc (A a ^{x)(t ) b ll c (x)) 



where 7 is defined by the horizon condition, [H], 

1 



Al 



" d v D, 



dN Al 2 
V2g 2 

dN A 
C A g 2 ■ 



2.11 



(2.12) 



(2.13) 



In noting this then it is elementary to reproduce Gribov's original gap equation by integrating 
the mixed propagator over the loop momentum p of ()2.9|) to give the MS result, P3IHI) 



1 



C A 




a + 0{a 2 ) 



(2.14) 



where a = g 2 /(16-K 2 ) and /x is the renormalization scale introduced to ensure that the coupling 
constant remains dimensionless in d dimensions when dimensional regularization is used. We 
have also absorbed the usual factor of 4-7re - ^ into fj? where 7 is the Euler-Mascheroni con- 
stant. Although the gap equation is divergent the infinities are absorbed in the renormalization 
constants defined by 



Wo 



Z A A a » , 

Z-^w , 9o 



= VZ c c a 
Z g 9 , 7o 



i ab = 
Z 7 7 



bf 



/-< 



ab 



(2.15) 
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where denotes a bare quantity and g and 7 are the running coupling constant and running 
Gribov mass respectively. We note that it was shown in [3H1 135] that in the Landau gauge 



Z c = Z, 

and, in our conventions, 

Explicitly at three loops one has in MS 

n 

7 7 (a) : 



(z A z c y 



Z„\/Za 



-1/4 



(2.16) 
(2.17) 



[16T F N f -35C A ]— + 



+ 



280C A T F N f - 449C1 + 192C F T F N f 

(486C(3) - 75607)Ci + (89008 - 15552C(3))C| T F N f 
+ (20736C(3) + 19920) C A C F T F N f - 12352C A T$N? 



a 

192 



3456C F T F N f - 844&C F T$N} 



6912 



0(a 4 



(2.18) 



Though the four loop MS expression for the colour group SU(N) is available in |45j . from 
knowledge of the Landau gauge gluon and Faddeev-Popov ghost anomalous dimensions. 



In the gap equation was used to show that in the infrared the ghost propagator is enhanced 
and behaves as jzzyz as p 2 — * 0. This was obtained by computing the single one loop correction 
to the ghost 2-point function and noting that in the p 2 expansion of the diagram the leading p 2 
correction was precisely the 0(a) part of (|2,14|) for all fi. More recently (|2,1|) has been used to 
compute the two loop correction to the horizon condition, (|2.12l) . and the ghost 2-point function. 
In [lUj the MS correction for massless quarks, is 



1 



C A 



+ 



"5 3. (C A1 ^ 

/2017 _ 11097 K 65 (C^\ 35_ / (C^ 

A \ 768 2048 256 sv ' 48 \ ^ J 128 \ y /i 4 

1137 n;,t n s 205tt 2 \ 
+ 7tt^v5C(2) 



2560 



512 



12 



7T 

+ -r 



+ 0(a 3 



(2.19) 



where S2 = (2v / 3/9)Cl2(27r/3), Ch{x) is the Clausen function and £(n) is the Riemann zeta 
function. Interestingly the 0(p 2 ) part of the 0(a 2 ) correction to the ghost 2-point function 
matched precisely with the 0(a 2 ) part of (|2.19j) for all \i so that in (|2.1jl the ghost propagator is 
still enhanced in the infrared limit. This is one of the conditions of the Kugo-Ojima confinement 
mechanism, |46j . which ensures that the Gribov-Zwanziger Lagrangian is an important tool for 
exploring the infrared structure of QCD from an analytic point of view. 

In the context of this Faddeev-Popov ghost enhancement at two loops in (|2.1j) we have 
also studied the structure of the 2-point function of the anticommuting Zwanziger ghost field 
ujf in the same infrared limit. Although this is a spin-1 field it has some similarity with 
the corresponding Faddeev-Popov ghost field c a since they both have the same wave function 
renormalization constant. It transpires that in evaluating the 1 one loop and 31 two loop 
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Feynman diagrams contributing to the uj^f 2-point function the same phenomenon occurs as 
the Faddeev- Popov ghost. The two loop gap equation, (|2.19|) . emerges in the same way in the 
p 2 expansion of the 2-point function so that the uj^ ghost is also enhanced in the infrared. In 
other words as p 2 — > the propagator behaves as rather than having the 4? behaviour 
of the initial propagator. It is worth remarking that as far as we are aware the discussion 
of |46j centred on a theory which involved only Faddeev-Popov ghosts and not the extension 
with the extra Zwanziger ghosts. Whilst in principle one would have to revisit that analysis to 
establish that the Kugo-Ojima mechanism still holds for (|2.1|) . the fact that the LO a ^ propagator 
has precisely the same structure of the Faddeev-Popov ghost would seem to imply that the 
Kugo-Ojima confinement criterion is still valid for (|2.1[) . Indeed it is worth observing that in 
the specific calculation, from the form of the u)^ propagator and its interaction with the gluon, 
the Lorentz structure of the lo^ 2-point function actually factors off and effectively leaves the 
Faddeev-Popov ghost 2-point function. 



3 Formal derivation of propagators. 

We now turn to the computation of the one loop corrections to the gluon propagator. However, 
given the nature of the problem we will also deduce information on the mixed propagator and <fi a ^ 
propagator simultaneously. First, we recall the normal procedure to determine loop corrections 
to a field where there is no 2-point mixing. Essentially one computes the corrections to the 2- 
point function and then inverts this expression truncating at the order in the coupling constant 
one is interested in. For (|2.1j) the procedure is the same. However, one is dealing with a 2 x 2 
matrix and has to first determine the corrections to each 2-point function which comprise the 
elements of the matrix. Once this has been achieved then the matrix can be inverted and all 
the propagators deduced from the appropriate elements. Several comments are worth making 
at this point. If one considers the location of the gluon propagator in the inverted matrix of 
2-point functions its structure is derived from two pieces. One is the determinant of the matrix 
which by definition involves all the elements of the original matrix and which will contribute to 
the denominator of the gluon propagator. The other piece will be located in the numerator of 
the final propagator correction and will also involve pieces from the Zwanziger ghost 2-point 
function. In other words the fields implementing the horizon condition will form a central, and 
as will be seen later, a crucial role in the behaviour of the gluon propagator and effective coupling 
constant in the infrared. It is worth recalling that these properties only result as a consequence 
of a non-zero 7 and hence the presence of a mixing term. 

We have described our strategy in a general formal way at the outset since the practicalities 
of the inversion is technical. First, one is dealing with spin-1 fields which require a gauge fixing 
term to have a non-singular inversion in momentum space. As this is standard we note that 
one can either introduce a gauge fixing term as we did in the previous section before setting 
a = to specify the Landau gauge we are interested in, or one can work in the Landau gauge 
directly and factor off the common transverse projector P lJiU (p) from the 2-point matrix. We 
choose to do the latter. However, another complication of the inversion lies in correctly taking 
account of the group theory structure of the 2-point function corrections. For the gluon itself 
and the mixed 2-point functions the one loop corrections have the same group structure as the 
Lagrangian term. However, as the 0^ 6 -field carries two group indices then its 2-point function 
will involve four colour index objects. If we denote these indices by the set {a, b, c, d} then we 
take the basis of independent rank four objects which can arise from the one loop diagrams as 

^ S aC 5^ 5 ac ^5^ c 5 a ^5 C( ^ J ace jbde jabe jcde ^abcd^ ^ -q 



s 



where the Jacobi identity excludes f ade f bce from being independent. Though we could equally 
have chosen either of the other two combinations. The object d\ bcd is defined as 



mbcd 

a a 



1. 



6 



-Tr [T1T^T C ,T 



A^A ± A ± A 



(3.2) 



where {T a ) bc = — if abc is the adjoint representation of the group generators and d°^ cd is the 
totally symmetric trace of four group generators in the R representation, 47 . To proceed with 
the inversion of the mixed 2-point function this colour structure has to be taken into account in 
the $> ab 2-point function and its corresponding inverse element. As this is a formal discussion we 
will not cloud the inversion by including the explicit results of the loop integrals at this stage 
but take the general structure of the matrices we are interested in as follows. The matrix of 
2-point functions at one loop is 



p 2 d ac 
jcab 



_~2 jacd 
_p2fiacfibd 



X5 a 



Uf 



acd 



4~ ^ jy^ jcab Q^ac^bd _|_ j^/ jace jbde _|_ jabe jede _j_ globed J ^ 4~ 0{tt ) 

(3.3) 

with respect to the basis j-^A^, <j!>™ 6 j where we have factored off the Lorentz structure. The 
quantities X, U, M, Q, W, R and S represent the one loop corrections. Similarly the Landau 
gauge inverse will be of the form 



+ 



\{p 2 ) 2 +C A ^\ 

( A5 cp 



cpq 



fped — \s cp 5 dq + 



7 fedr fpqr 

" ' r p 2 [{p 2 ) 2 +C a ^] j 1 



Cf 



cpq 



\ Ef pcd G5 Cp 5 dq + Jf c 'P £ f d ^ + K fcdepqe + L( fdl 



cdpq )a + 0{a 2 



(3.4) 



where we have included the propagators from the previous section and the quantities A, C, E, 
G, J, K and L will depend on the one loop corrections defined in (|3.3|) . Given these forms it is 
straightforward to check that 



.4 
C 
E 
G 
K 



1 



[(p2)2 + Ca7 4]2 



(p 2 ) 2 X - CaiYU - CaTP M + Ca1 A (Q + C A R + \C A W) 



[(p 2 ) 2 + c A7 4 ] 2 
i 

WY + Ca^ 
Q 



j 



l 



[(p 2 ) 2 + C A7 4 ] 2 

„4 



7 V X - C A <fM + (p 2 ) 2 U - tV (q + c a r + \C A W) 
jVX - C A J 4 U + {j?fM - 7 V (Q + C A R + \C A W) 

w _ s 

(p 2 ) 2 ' (p 2 ) 2 
7 4 X + 7 2 p 2 U + tVM + ( P 2 ) 2 R - 7 4 (Q + §CW) 



+ 



7 



(3.5) 



by ensuring that the O(a) term of the product of (|3.3j) and 1)3.4(1 vanishes. As a check on this 
inversion we note that C = E if M = f7. Also in determining the 0(a) correction to the mixed 
propagator of (|2.7j) C and £7 have to be divided by v2 as was noted after (|2.8|) . We have derived 
these corrections at the formal level to indicate how involved the one loop corrections to the 
propagators are. However, given the intricate relation with the 2-point function corrections it 
will transpire that they will only be of use at this formal level. 
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4 One loop corrections. 



We now detail the computation of the one loop integrals. First, the required Feynman diagrams 
are generated with the use of the Qgraf package, [3H!- For the gluon and 4> a ^ 2-point functions 
there are eight and two diagrams respectively whilst there are two diagrams contributing to 
either of the mixed 2-point functions giving a total of fourteen. Although we are primarily 
interested in examining the p 2 — ► limit of the propagators to ascertain whether the one loop 
gluon propagator still vanishes as it does at one loop, we have decided to evaluate the one loop 
correction exactly as a function of p 2 . This requires the explicit evaluation of the one type of 
master integral which is of the form 



Ix{p,m\,rnl\a,ff) = \ 

Jk 



k [k 2 + m\ ) a [(k - p) 2 + ml}/ 3 



(4.1) 



where at one loop a and (5 are integers and the mass arguments take any combination of values 
in the set {0, i^J CaJ 2 , —WC-aI 2 }- This form, (|4.1|) . derives from rewriting the propagators with 
the Gribov structure p 2 /[(p 2 ) + Ca7 4 ] using partial fractions such as 



[(p 2 ) 2 + CU 7 4 ] ' 2 I [p 2 + ijC~ A ~i 2 \ 




(4.2) 



producing standard propagators but with a width. However, it is worth noting that in any 
final expression we derive, the answer must be real given that the initial integrals are real 
which provides an internal check on the computation. On a practical note we record that the 
decomposition (|4.2j) is readily implemented in a computer algebra programme written in the 
symbolic manipulation language Form, [1^1 • Indeed the complete calculation we describe here 
has been performed automatically with Form. 

Whilst the type of master integral, Ii(p,m 2 ,m 2 ; 1, 1) has been studied and exploited many 
times we note that the key difference here is the presence of the complex mass. However, in the 
exact evaluation of our integrals we note that we use the formal results for Ji(p, m\ , m 2 ; 1, 1) 
with real mf before analytically continuing to the values ±iy / Cvl7 2 we are interested in when 
the masses are non-zero. For completeness we note that the results for the two central integrals 
we use are, expanded to the finite parts, 



Ii(p,iVc^j 2 ,iVCA~l 2 ]l,l) 



1_ / 16g A7 4 
2 n ^ (p 2 ) 2 




x tan 



V2 



-1/2-1 
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x tan 



(p 2 ) 2 
V2 



+ 1 



1/2 



1 + 



+ I 



16CU7 4 



2^2 



16Ca7 4 

(p2)2 

1/2 i ta 

2 



-1/2- 



( WC A ^ 
{ (P 2 ) 2 , 



2 
In 



In 



i + W(i + 



16CA7 4 



,2\2 



1 + 1 1 + 



16C A 7 4 



(p 2 ) 2 



1/2 



+ O(e) 
(4.3) 



and its conjugate, and 

h( Pl iVC^7 2 ,-iVC^l 2 -AA) 



- + 2 

e 



[p 2 - 2^^172] 1 / %\[Cjci 



2p 2 



In 



b 2 + 2i^ 7 2 ] 
2p 2 

\/[4Ca7 4 - (p 2 ) 2 



In 



-i%/Cl7 2 



tan 



P 



V[4C A7 4 - (p 2 ) 2 



P 



+ 0(e) . 



(4.4) 



When one of the masses mf is zero, we use the result 



k k 2 [(k - p) 2 + iVU^'j 2 



+ 2 



In 



i^Cll 2 



[p 2 + iVClj 2 } 



I 1 



In 



\p 2 + i^fCl-j 2 
iJCll 2 



+ 0(e) (4.5) 



and its conjugate. Finally, for completeness the one loop vacuum bubble is 



Ik 



k [k 2 + i^Cli 2 ^ 



r(i - ¥) 

(47r)d/2 



(^7 2 )" 



a!— a 



(4.6) 



which, together with its conjugate, is expanded in powers of e and the resulting logarithms 
treated with 



~1*(^7 4 ) + £ 



(4.7) 



and its conjugate. We now discuss the 2-point functions themselves but first note that they have 
to be renormalized. These infinities are absorbed by the wave function renormalization constants 
given in section 2, and we follow the procedure used in !")(.). to remove them in an automatic 
calculation. Though we note that the mixed 2-point function is in fact finite as a consequence 
of the Slavnov- Taylor identity relating Z 7 to Za and Z c . Given the rather involved form for 
the master integrals, (|4.Hjl and (|4.4j) . the full explicit expressions are large and do not serve 
to illustrate any major points which can be accessed by other methods we will discuss. Since 
we are primarily interested in the zero momentum limit, we will concentrate on extracting that 
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behaviour. There are two ways of doing this. One is to expand the explicit expressions in powers 
of p 2 and truncate at the appropriate term. A second way is to return to the Feynman diagrams 
themselves and carry out the expansion of each integral in powers of p 2 . This is known as the 
vacuum bubble expansion and is one method to renormalize a quantum field theory. Though of 
course applying this technique here would require the term which is finite in e. We have chosen 
to do this calculation in addition to the exact calculation for two reasons. First, it provides an 
alternative check on the explicit result where the expansion of the exact result in powers of p 2 
must agree with the vacuum bubble expansion. Second, if one were to go beyond one loop to 
carry out a two loop analysis, then the exact evaluation of the two loop integrals with masses 



in the set {0, i\JCai 



'CaI } would be a formidable task even if all the massive master 



integrals were known in closed form explicitly before analytic continuation. The vacuum bubble 
expansion provides a more practical and efficient approach. To proceed with the expansion, the 
propagators of the master integral Ii(p,m 2 ,m 2 ;a, (3), for example, can be expanded in powers 
of p 2 where the p-dependent propagator is expanded recursively with 



p) 2 + m 2 ] 



+ 



2kp — p 2 



[k 2 + to. 2 ] [k 2 + to 2 ] [(k — p) 2 + to 2 ] 



(4. 



which is readily implemented in FORM. The truncation criterion is that the 2-point functions 
themselves are O ((p 2 ) 2 ). Given the fact that (|2,lj) is renormalizable then the O ((p 2 ) 2 ) terms 
will be e-finite. However, in performing this bubble expansion we cannot apply the identity 
(|4.8|) when all the masses in a Feynman integral are zero. This is because such an integral has 
its momentum dependence predetermined by (p 2 ) r where T is the dimension of the integral. 
As the momentum dependence is fixed due to the masslessness of the graph one needs to first 
isolate such integrals prior to applying the vacuum bubble expansion to the remaining terms in 
the decomposition of the overall Feynman graph. Once the vacuum bubble expansion has been 
performed it is elementary to replace the vacuum bubbles with the general result 1)4. 6j) . 

Given this algorithm it is straightforward to implement it in Form and we record the p 2 
expansion of each of the 2-point functions. We find 



(A*J-p)A>(p)) 



(A«(-p)& c (p)) 



cab 



P 2 + 



+ 0((p 



,2\2 



eabc 



1 + 




§ ac 6 bd 



64 

+ O ((P 2 ) 2 ) 



jabcd P 

32 C A 



(4.9) 



which determine X, M, U, Q, W, R and S. We note that the expansion using Q4.8|) agrees with 
the expansion of the explicit expressions on the right side of (|4.3J) . (|4.4[) and 1)4.6(1 . From (|4.9|) 
we can deduce several interesting properties of the p 2 — > limit. First, it is clear that at one 
loop the gluon propagator, (|3.5|) . vanishes as the momentum vanishes. This is because when one 
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substitutes the explicit values for X, M, U, Q, W, R and S from (|4.9|) into the expression for 
A in ()3.5|) then one finds 



p2 


^ In 

8 


fc A1 A \ 


215~ 


Ca1 a 


K ^ ) 


384 



a + 0[(p 



,2\2 



+ 0(a 2 



(4.10) 



The leading terms in this expression derive from only Q, R and W. This vanishing of the gluon 
propagator at one loop is then consistent with general expectations and lattice results. See, for 
instance, |51|. Further, if one examines the 8 ac 8 bd channel of the cj)^ propagator the one loop 
mass gap condition (|2.14j) alters the p 2 dependence in a way similar to what occurs in both the 
c a and uff 3 ghost propagators. Moreover, this is the channel which appears in the Lagrangian 
itself and this implies that if there was no complication from the mixing then there would be an 
enhancement in that channel. 

Finally, we note that the full explicit expressions for the one loop corrections to the 2-point 
functions are recorded in appendix A. Whilst these can be substituted into the formal expressions 
for the propagators, (|3.5|) . the explicit expressions do not lead to any further insight. Though 
the formal functional form may be of use in obtaining better parametrizations of lattice and 
DSE data in the low p 2 region. However, for the remaining fields which are the Faddeev-Popov 
ghost and the anticommuting Zwanziger ghost, the propagators are 
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where the ghost form factor is defined by 

(c a (p)c\-p)) 
Using the gap equation, (|2.14|) . we have 
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C A 1 A 
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2\2 
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(4.11) 



(4.12) 
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(4.13) 



So that D- l (p 2 ) is 0{(p 2 ) 2 ) as p 2 



0. For uj® b it turns out that defining 



«(pK»(-p)) 



cd I 



gac £>d {P 2 



(4.14) 



then 



DM 2 



D c (p 2 



(4.15) 
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to one loop and no longitudinal component is generated. So the lo^ ghost enhancement actually 
follows trivially from the equivalence with the Faddeev-Popov ghost for all momenta. This 
follows as a result of our earlier observation of the factoring of the Lorentz structure from the 
2-point function. Given there is a similar enhancement at two loops from the vacuum bubble 
expansion, it is tempting to speculate that not only does 1)4.15(1 hold at two loops but maybe 
also to all orders. Whilst we have noted that there is a similar type of enhancement of Q in 
the (^f-propagator, it is also turns out that Q = — (D c (p 2 )) = — (D u (p 2 )) at one loop. 
Finally for the quark propagator, where the quarks are massless, defining the form factor by 



(ip ll (p)^ jJ (-p)) = i5 l3 5 



u 



(4.16) 



then 
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We note that in the limit p 2 — ► then 
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C F a + 0(a 2 



(4.17) 



(4.18) 



so that the form factor tends to a constant at zero momentum. 



5 as freeze-out. 

Having examined the zero momentum limit of the 2-point functions and form factors, we now 
apply the results to the problem of the value of the effective strong coupling constant in the 
same limit. In both lattice and DSE analyses one can study the value of a renormalization 
group invariant effective coupling constant, denoted by alf(p 2 ), which is defined as a result of 
the renormalization properties of the gluon ghost vertex in the Landau gauge. In terms of the 
gluon and Faddeev-Popov form factors, Da(p 2 ) and D c (p 2 ), where 

(A°(p)AU- P )) = 5 ab ^lp^p) (5.1) 

the effective coupling is 

a|V) = a^)D A {p 2 )(D c {p 2 )f (5.2) 

where a(/j,) is the running strong coupling constant and a = g 2 /(4-7r). Due to the gluon ghost 
Slavnov- Taylor identity, |52j . there is no contribution from the vertex form factor. Thus, given 
that we have one loop expressions for Da(p 2 ) and D c (p 2 ) in the p 2 — > limit, it is straightforward 
to examine the structure of 1)5.2(1 in the same limit. 

First, if we consider the expression with the tree level form factors 

D A(P 2 ) = p#P^4] + °(«) > D ^P 2 ) = 1 + °( fl ) ( 5 - 3 ) 
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then clearly 



«!V) 



a((j,)(p 



2\2 



[(p2)2 + Ca7 4 



(5.4) 



which vanishes as p 2 — > 0. However, if one includes the one loop corrections then with 



D C ( P 2 



1 + 



5 3. C A7 A 
8"8 ln h^ 
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+ 0{p 4 



C A a + 0{a 2 ) 
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(5.5) 

where we have expanded the propagator to the next term in the momentum expansion, we have 
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1 + C A ( | In 



CA7 4 (M) 



5 j 7rp^ 

8 "T" Sy^T 2 ^) 



a(/x) 



(5.6) 



In this expression we have expanded the denominator factor of D A (p 2 ) and included it in the 
0(p 2 ) part of the numerator. Though in fact it will not contribute to the final value of a|?(0) 
which is derived from the displayed numerator and denominator factors. Before taking the limit 
one can use Gribov's gap equation, (|2.14j) . to enforce the ghost enhancement in the ghost form 
factor which leaves the denominator as 0(a 2 ). This coupling constant dependence in fact cancels 
with similar factors in the numerator. One comes from the a(fj,) in the definition and another 
one derives from the numerator factor when the gap equation is used. Unlike the cancellation 
in the ghost propagator to give ghost enhancement, in the numerator the numerical term at 
one loop does not match that of the gap equation and one is left with an 0((p 2 ) 2 ) term in 
the numerator. This momentum dependence then cancels the (p 2 ) 2 term on the denominator. 
Finally, the running Gribov mass parameter, 7(/u), also cancels with one factor of 7 deriving 
from the determinant of the 2-point functions and one from the p 2 / 7 2 term of the expansion of the 
Faddeev-Popov ghost 2-point function. Since the running coupling constant and dimensionful 
quantities also cancel one is left with a pure number when p 2 — > 0. We find 



af(0) 



50 



3irC A 

Specifying various gauge groups of interest, we have 



af(0) 



5(7(3) 



and 



af(0) 



SU(2) 



1.7684 



2.6526 



(5.7) 



(5.8) 



(5.9) 



This is a novel property of the Gribov- Zwanziger Lagrangian. From studying the calculation 
in the way we have presented it, we note that the gap equation was essential in removing the 
tree terms of the form factors, as well as providing factors of a(n) and 7(/x) which precisely 
cancel. Indeed the role of 7(/x) seems crucial since it provides the balancing dimensionality for 
the momentum which eventually vanishes. More significantly it is interesting to note that the 
cancellation between the numerator and denominator factors is between pieces which involve the 
two types of ghost. In other words the anticommuting Faddeev-Popov ghost and the commuting 
Zwanziger ghost, (j)^ . The contribution from the former, which is required for the usual gauge 
fixing, is clearly demanded by the definition but the latter arises through the implementation of 
the Gribov horizon which is, of course, an infrared contribution originating from the ambiguity 
of the gauge fixing procedure. 

Concerning the actual numerical value of a|?(0) in SU(3), (]5.8[) . it is worth contrasting with 
estimates from other methods. From a survey of articles (which is by no means exhaustive) 



15 



it appears that the non-zero values fall into three classes. These are when a| F(0) is deduced 
from experimental data (P), from numerical studies such as lattice or DSE and from more 
analytic approaches, (A). Currently there is not a consensus of values from these three divisions 
which are summarised in Table 1 where the result of |TS] is the central value of the range the 
authors quote. For instance, the DSE value differs from the phenomenology based approach by 
an order of magnitude; with analytic approaches, including this article, bridging between the 
two extremes. However, in this context we note that 1)5.7(1 should be viewed as a qualitative 
result rather than quantitative. First, the computation, whilst simple in its derivation, is only 
at one loop. Clearly there are higher order corrections to the form factors. Indeed it is not 
immediately apparent how the 0(a 2 ) corrections would conspire to cancel to leave a result 
independent of a(/z) and as was the case at leading order. However, it is also not clear 
whether comparison of a|?(0) from different methods has any meaning in the first instance since 
it corresponds to the evaluation of a quantity defined in perturbation theory but evaluated in 
the infrared limit. Though the lattice and DSE estimates are for the same quantity as (|5.2j) 
in the Landau gauge. Moreover, 1)5,2)) depends on the ghost form factor which would not be 
an object immediately accessible experimentally. In this context it is worth commenting on a 
recent proposal of |321 where an effective coupling constant is defined from other vertices 
in the Gribov-Zwanziger Lagrangian. For instance, it was argued that one could examine the 
form factor of the triple gluon vertex, Daaa(p 2 ), arid derive an effective coupling with Da(p 2 ), 
defined with respect to some momentum configuration for the external gluons, with 

afW/) = a(n)D AAA (p 2 )(D A {p 2 )f (5.10) 

where D AAA (p 2 ) is necessary due to the lack of a Slavnov- Taylor identity comparable to that 
of the gluon ghost vertex. To have a freezing of ot^ AAA (p 2 ), given the behaviour of D A (p 2 ) 
established here, then D AAA (jp 2 ) would have to diverge in such a way that the momentum 
dependence cancels in the zero momentum limit, [32( I33j . Whilst this is essentially a dimensional 
observation, |32l 133 j . it does not necessarily follow that the same numerical value for the frozen 
effective coupling would emerge. Given that we now have a renormalizable localized Lagrangian 
implementing the Gribov horizon, one natural question which could be considered is whether 
this behaviour for the triple gluon vertex could be determined in the p 2 —* limit. Although 
such a calculation is beyond the scope of the current article, we note that on renormalizability 
grounds D AAA (p 2 ) would have to be finite as p 2 — > 0. Therefore, if one is to extract a value for 
a !s^L4A(0)> via t ne suggestion of [321 ES], it would seem to us that the behaviour is truly driven 
by some non-perturbative mechanism. 
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Table 1. Non-zero estimates of a|F(0). 



Although we have given a summary of the non-zero finite values for a| (0) it should be noted 
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that various lattice and DSE studies find a value of zero as p 2 — ► 0, (see, for example, [SUEZ!), 
for the effective coupling defined from the gluon and ghost form factors. Further, lattice and 
DSE studies of effective couplings defined from other vertices also find a zero value in this limit. 
For instance, see |35U36| . Therefore, it would appear that currently there is no common view of 
the precise zero momentum behaviour. However, these latter studies are of the effective coupling 
defined from the triple gluon and quark gluon vertices respectively. So if our argument about 
the finiteness of Daaa(p 2 ) as p 2 — ► from 1)2 .1JI is valid, then the vanishing of Da(0) would 
imply that &saaa(®) = ^ which would appear to be consistent with the analyses of this vertex, 



The situation for SU(2) is similar, though of course the only results available are from lattice 
and DSE studies. Briefly, in the same way that certain SU(3) lattice and DSE computations 
give similar non-zero values for a|?(0) the same is true for SU(2) with the common freeze-out 
value of a| ff (0) = 5(1), [En], or 5.2, 51 . Whilst this is significantly larger than our value of 
2.65 there is at least the consistent observation that our SU(2) value is larger than the SU(3) 
value. Though from the explicit expression this is primarily due to the denominator color factor 
in 1)5.71) . Taking the ratio of the SU(2) and SU(3) values for a|?(0) from the DSE analyses one 
finds the ratio 1.75 in comparison with the ratio of 1.5 of (|5.7]> . 

Finally, in comparing IJ5.7JI with DSE estimates it is worth noting that the latter are invariably 
computed in MOM schemes which are mass dependent renormalization schemes. By contrast 
(|5.7|) has been deduced in the MS scheme which is a mass independent scheme. Whilst ultimately 
the definition of a^(p 2 ) is a renormalization group invariant, if the quantity (|5.2[) has a meaning 
in the infrared then the full value will be scheme independent. Therefore, computing the two 
loop corrections to Q5.7JI would be useful to see if there is convergence to a higher value. In 
discussing some of the issues concerning the freezing from a theoretical point of view, it is worth 
noting the one loop study of :55^ where the infrared behaviour of an effective coupling defined 
via the quark gluon vertex with massive quarks in an arbitrary covariant gauge in the MOM 
scheme. There the presence of the gauge parameter led to differing zero momentum behaviours 
indicative of the subtleties associated with analysing the infrared behaviour of a quantity which 
has only justifiable meaning in the ultraviolet regime. 

6 Power corrections. 

Having studied the effective coupling in the limit of zero momentum, using both the bubble 
expansion and the exact expressions for the 2-point functions, we can now examine the propa- 
gators and the couplings in another limit. In |431 144j . a numerical fit of the lattice computation 
of the effective coupling constant suggested that in a certain momentum range, the coupling 
deviated from the expected (perturbative) behaviour by a piece which could be parametrized by 
power corrections. Intriguingly these corrections were of the form 0(l/p 2 ) and not the expected 
0(1/ (p 2 ) 2 ) behaviour. In other words, [4~3l 14"4"| . 



An 0(1/ (p 2 ) 2 ) correction is motivated by the fact that to match the dimensionality of the mo- 
mentum dependence one needs a dimension four quantity. Given the fact that the operator 
product expansion and sum rule studies of suggest the condensation of the gauge invari- 
ant operator G < / lu G afML ', then there is an a priori clear candidate for the numerator which is 
{G® G ^}. However, in order to have an 0(l/p 2 ) correction, one would have to have a lower 
dimensional operator condensing. Since one such operator is ^A^A "^, this has led to the sug- 




(6.1) 
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gestion that this operator condenses, [15) . and consequently estimates for (^A^A a ^) have been 
derived from the operator product expansion, [S7j, and the local composite operator method, 
[5%] . Indeed it was observed in [23 EHJ that the perturbative QCD vacuum was unstable and the 
vacuum expectation value of such operators, and its non-local gauge invariant generalization, 
would play a significant role in understanding the true vacuum. Also, the condensation of such 
a dimension two operator was already discussed in the Coulomb gauge in |6*T] . Clearly \A a p .A a ^ 
is a gauge variant operator but the appearance of its vacuum expectation value in the fits of 
(431 144j is neither unexpected nor inconsistent with this since the effective coupling constant is 
gauge dependent in mass dependent renormalization schemes such as the MOM scheme used in 
|43| I44j . As a consequence of these observations concerning the apparent existence of a non-zero 
value for (^A^A afl ) there has been renewed interest in trying to understand the dynamical gen- 
eration of a gluon mass. However, given that the Gribov mass, 7 2 , can effectively reproduce a 
gluon mass, it is worth investigating whether one can mimic the power correction behaviour of 
|43| 144] from the Gribov- Zwanziger Lagrangian, (|2.1|) . Therefore, we have examined the integrals 
used to construct the form factor corrections exactly and expanded them in powers of y / C^7 2 . 



Though we remark that in the expansion in which we work yXJ 
rewrite one of the parts of an exact integral through the replacement 
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Therefore expanding the 2-point functions in powers of 7 2 , produces 
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and 
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-1 



(6.4) 



where we have neglected the term beyond the first 7 2 /p 2 correction in each case. For the gluon 
and Faddeev-Popov ghost the finite 7-independent pieces both agree with the result of the 
massless calculation. 

Interestingly the gluon propagator correction is 0(7 2 ) and not 0(7 4 ) as might have been 
suggested from the original gluon propagator which is a function of 7 4 . Although the correction 
is from the one loop term and involves a, it qualitatively introduces an effective mass for the 
gluon in this next to high energy limit. One comment concerning the choice of the sign of 7 2 
is worth making. Eliminating the auxiliary scalar field (f)^ from (J2.1|) means that Q2.1[) is an 
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even function of 7 2 . One could fix the sign of 7 2 by arguing that it has to produce a non- 
tachyonic effective propagator of the usual form from the binomial expansion. However, one 
could equally choose the mass to be tachyonic without upsetting the two loop gap equation and 
the ghost enhancement. Although this may not appear to be a physically sensible choice, we 
remark that the study of 62 suggested that introducing an effective tachyonic gluon mass into 
current correlators could reproduce experimental data more accurately compared to the case of 
no gluon mass. Whilst there did not appear to be any justification for such a tachyonic gluon in 
|62| . the freedom of choosing the sign of 7 2 , whilst still retaining the infrared properties already 
discussed, does appear appealing. 



In addition we can now examine the effective coupling constant, 1)5. 2[) . in the same limit. 
Retaining only the 0(7 2 ) piece of A from the 2-point functions, then 

2 3C a kV(W ( 1 \ 

D A(P ) = 1 g ^2~ a + ° [j^2 ) ■ ( 6 ' 5 ) 

In expressing the power corrections in this way, we are treating it as an expansion in powers of 
VCaJ 2 /p 2 since the appropriate factor of Ca is associated with the Gribov mass. Thus from 
the definition, and regarding the 7 independent piece as corresponding to a^ ert (p 2 ), then to one 
loop we have 

eff/ 2\ pert/ 2\ -, 9C A ^ VCaJ 2 „ / 1 \ (R R , 

«S (P ) = «S (P ) 1 g ^2~ a + [j^2 ) ( 6 - 6 ) 

in MS. From the 0(7 2 ) corrections one formally produces an 0(l/p 2 ) term as the leading 
power correction as opposed to an 0(l/(p 2 ) 2 ) one. Also the coefficient has the opposite sign 
to that of (431 144j and would be in keeping with our suggested choice of changing the sign of 
7 2 . However, since that calculation was in the MOM scheme in contrast to the MS scheme 
one loop computation here, this would not necessarily be sufficient justification for altering the 
sign of 7 2 . Again we merely regard it as an interesting observation. Further, we emphasise 
that, as in earlier sections, these are qualitative observations of this effective coupling constant 
deriving from a one loop calculation with the motivation being to appreciate the implications 
of the Gribov parameter in comparison with non-perturbative numerical analyses. Clearly such 
a power correction could equally well be produced by the presence of (^A'^A afl ) instead of the 
condensate (|2.12|) . Indeed there have been more recent studies, [23 ESI EU , of the effect a non- 
zero {\A^ l A a ^ 1 ) has in the Gribov- Zwanziger Lagrangian in the Landau gauge. We also remark 
that from the renormahzation structure, |38[ I39| . both \ A a ll A a ^ and 7 2 have the same anomalous 
dimensions. So that a combination of both these dimension two operator vacuum expectation 
values could be responsible for the 0(l/p 2 ) correction in the data of |43[ 144 j . At a more formal 
level it is worth commenting on the structure of the 0(7 2 ) correction in (|6,6[) . This effective 
coupling constant is by construction a renormahzation group invariant quantity due to the 
underlying Slavnov- Taylor identity for this vertex, [HE]- For the Gribov-Zwanziger Lagrangian 
this identity remains valid which can be observed, for instance, from explicit loop calculations 
of the relevant anomalous dimensions using (|2,1|) which reflect the ultraviolet position. As is 
evident from the freezing calculation, the numerical value of a|?(0) is clearly a renormahzation 
group invariant. However, it is not immediately apparent if the power correction, which is 
manifest in the remaining energy range, also retains this feature. To fully examine this situation 
is not as straightforward as it would simply appear. This is primarily due to the fact that 
the full renormahzation group for (|2.1j) has not yet been constructed. Ordinarily this involves 
parameters and quantities which are independent. By contrast, in (|2.1jl the Gribov mass, being a 
parameter of the theory which appears in the (ultraviolet) renormahzation group equation, is not 
independent by virtue of the Gribov gap equation. Therefore, to fully access the renormahzation 
group properties of quantities deduced from (|2.1|) it would seem necessary in the first instance 
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to develop the renormalization group for (|2.1j) subject to the Gribov gap equation constraint. 
This is beyond the scope of this article. 

If we also consider the quark propagator, then in the same limit we find 



By the same token that we argued that the power corrections appeared to generate a mass for 
the gluons one might regard this correction as a generated quark mass. Although the sign is 
opposite to what one observes in the gluon case but with magnitude 37rC^a/4, we emphasise that 
this is in the situation where there is no initial quark mass and therefore it does not correspond 
to the part of a fermionic propagator which ordinarily defines the mass term. 

Finally, if one accepts that there can be an effect from the Gribov parameter in studying 
power corrections in gauge variant quantities, then considering higher power corrections leads 
to the conclusion that one could have contributions from not only (G®G a,iV ) but also 7 4 in the 
0(l/(p 2 ) 2 ) correction. Therefore, in extracting estimates for (Gr" G ativ ) in such calculations, 
one would in principle need to make allowance for the potential presence of additional 7 4 type 
terms. For the computation of gauge invariant quantities, there ought not to be any contribution 
from the Gribov mass. 

7 Discussion. 

In examining the one loop corrections to the propagators in the Gribov- Zwanziger Lagrangian 
several interesting features have emerged. First, we have shown that the gluon propagator 
vanishes in the p 2 — > limit which is consistent with DSE and some lattice studies. More 
significantly we have seen to what extent the gap equation satisfied by the Gribov parameter 
underlies the infrared structure of the propagators and effective coupling constant in ()2.1[) , In 
it was used to ensure that the ghost propagator was enhanced in the infrared. However, 
here it has been used to show that the u;° b ghost propagator is also enhanced at two loops, 
and that the behaviour of the tree channel of the </>^ b propagator is also altered in the infrared 
at one loop. More significantly the gap equation played a central role in the freezing of the 
renormalization group invariant effective coupling defined via the gluon Faddeev-Popov ghost 
vertex and has given some insight into how a non-zero finite value emerges. For instance, in the 
definition (|5.2|) if the ghost propagator does not satisfy the gap equation then the factors of p 2 
in the zero momentum limit will not match between the numerator and denominator to leave a 
momentum independent quantity. Any remaining momentum dependence would dominate the 
p 2 — > behaviour. Since there appears to be a discrepancy in various numerical approaches 
between a non-zero and zero value for a|?(0), it may be due to the Kugo-Ojima condition not 
being fully and precisely satisfied. Although our one loop calculation has led to a freeze-out value 
which is different from that from DSE and certain lattice studies, and is regarded as qualitative 
in the sense described earlier, it would be interesting to see what effect the two loop corrections 
have on the numerical estimate. At one loop the result is clearly independent of the number of 
quarks but quark loops will be present at two loops. Recently in [EE] a comparison was made 
between a calculation of (|5.2|) in quenched QCD and one where dynamical fermions were used. 
Although both computations suggested a| ff (0) = 0, the peak of the plot of a^(p 2 ) is lower in 
the dynamical case. Though we note that intriguingly the peak in the quenched case appears to 
be around 1.7 which is similar to ()5.8|) . In concluding this aspect of our analysis it would be fair 
to make the general observation that the full resolution as to whether the quantity (|5.2|) freezes 




(6.7) 
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at a zero or non-zero value in the infrared has not yet been achieved and our value is merely 
another contribution to that debate. 



Another aspect of this study has been to appreciate the role played by the extra Zwanziger 
ghosts. There are various ways of viewing them. In one both and lo°^ can regarded as 
being on the same footing from the way they were originally introduced to localise the Gribov 
problem. However, from the equation of motion for <#J° ', which was used to implement the 
horizon condition, ^ is clearly not unrelated to the original gluon field itself and in some sense 
could be regarded as part of a more general spin-1 field such as 

Al = A" + ^-^-A^ (7.1) 

where A is a constant, and which could be regarded as taking account of some of the more global 
aspects of the gauge fixing. Indeed the freezing of the effective coupling was dependent on the 
value of the ^ 2-point function. So in some respects the infrared behaviour is guided by the 
inherent non-locality of a more general spin-1 field. Whilst this field A^ is clearly non-local it 
is similar to the second term in the expansion of the gauge invariant spin-1 field of (HZj which 
is also non-local. In _S7| the non-local structure resulting from the gauge invariant mass term 
in the Lagrangian was used to construct a vortex solution which underpinned the confinement 
process. Given that there is now numerical evidence for vortices present on the Gribov boundary 
(in the Coulomb gauge [Ml ) , it might be possible to construct an explicit vortex solution in the 
Gribov-Zwanziger Lagrangian and study their dynamics in relation to the role they have in 
confinement. In contrast to the non-local mass term of |)7], which is non-localizable since it 
would require an infinite number of extra fields, the non-locality of the Gribov formulation of 
non-abelian gauge fixing is localizable since only the finite set of fields {(j>f, 4> a J>} and {ujf,ijf} 
are required. A heuristic way of viewing this is to formally eliminate (\>°^ from (|2.1|) . and ignoring 
all anticommuting ghost fields for the moment, then one has 

rGZ _ 1 n a r<aiiv C A7 4 A g 1 Aa a dN A ~/ 4 

- 4 " — A »~&W A + w 

which involves only one non-local term. In momentum space the quadratic part of (|7.2|) clearly 
leads to the Gribov gluon propagator. Returning to A^, one can now view the anticommuting 
Faddeev-Popov and Zwanziger ghosts from a different point of view to their original role, if one 
ignores their spins. In other words it seems that one can associate a ghost field with each term 
of A^ and, as we have shown here, each of these ghosts have similar infrared properties. For 
instance, both their propagators and form factors are equal at one loop and enhance at two 
loops in the p 2 — > limit. Also the renormalization group invariant couplings defined from their 
interactions with the gluons both freeze at zero momentum to the same value. In this respect, 
aside from their differing spins they could be regarded as the first few terms of a ghost multiplet 
where the subsequent terms would be associated with higher order terms of A° which would 
involve higher dimension non-local operators. Such operators might be the ones necessary to 
give the potential freezing in the effective coupling constants derived from the triple and quartic 
gluon vertices as suggested in [321 I33| or the restriction to the fundamental modular region. 

Finally, we note that we have concentrated throughout in this article on studying the Gribov- 
Zwanziger Lagrangian in the Landau gauge only. However, it would be interesting to repeat 
the present analysis for other covariant gauges such as the other linear covariant Lorentz gauges 
and the maximal abelian gauge. Indeed for the former gauges there has been preliminary work 
in this direction, |H9[ I7()j. and it would be interesting to see, for instance, whether the power 
correction behaviour changed significantly and if the variation in magnitude of the residue of 
the 0(l/p 2 ) correction, if any, could be measured, say, on the lattice in a variety of gauges. 
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Moreover, in the maximal abelian gauge it is possible to define an effective coupling constant 
from the Slavnov- Taylor structure of one of the gluon ghost vertices, [70]. Therefore, if one 
obtained the same freeze-out value for this effective coupling, it could be evidence for the abelian 
dominance hypothesis. We close by remarking that our study implies that the Gribov-Zwanziger 
Lagrangian has established infrared features which are generally not inconsistent with other 
methods. Specifically, these involve quantities such as ot e ^{p 2 ) which do not ultimately involve 
the running parameters g((i) and y(fJ>)- Therefore, it would be interesting to apply (|2.1|) to other 
(infrared related) problems where such dependence on the running parameters cancels. 
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A Explicit 2-point functions. 



In this appendix we record the explicit forms for the one loop corrections to each of the 2-point 
functions. Using the compact notation of (|3.3|) we have 
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We note that in these expressions the explicit [i dependence appears only in a small set of terms, 
which can clearly be identified with the piece of the 2-point function deriving from the ultraviolet 
contribution. 
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